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ON THE LAGKANGTAN METHOD FOR 
STEADY AND UNSTEADY FLOW 


W. //. Hui and II. J. Van lloeatel* 

Department of Applied Mathematics 
University of Waterloo 
Waterloo, Ontario, Canada. 


ADSTIiACr 

A new and general Lagrangian formulation of fluid motion is given 
in which the independent variables are three material functions and a 
Lagrangian time.,, which differs for different fluid particles and is distinct 
from the Eulcrian time. For steady flow it requires only three indepen* 
dent variables - the Lagrangian time and two stream functions • in con* 
trast with the conventional Lagrangian formulation which apparently 
still requires .four independent variables for describing a steady flow. 
This places the Lagrangiau formulation for steady flow on the same foot* 
ing as the Eulcrian. For unsteady flow, the new formulation includes the 
conventional formulation as a special case when the Lagrangian time is 
identified with the Eulerian time and when the material functions are 
taken to be the fluid particle’s position at some given time. The distinc* 
lion between the Lagrangian and Eulcrian time, however, is found useful 
in applications, e.g. to problems involving a free boundary. 


'Pitmit addreia: Depaitracbt of Applied Mitkenutio, Uiinmljr of Weitem Ontario, 
Londoa, Oatario, Canada. 


1. INTRODUCTION 


It it well-known that there exist two general methods of description of fluid 
motion. The Eulerian method aims at describing the flow variables, e.g. velocity v, as 
functions of the position x ■ (*,V,z) of a fluid particle at time t, vis., v * v(z,y ,*,(). 
In contrast, the Lagrangian method aims at describing the motion history of each fluid 
particle whose itame is (<i,t,e), vis., x ■ x(o,k,c,t). Conventionally, the Lagrangian 
variables (a,6,e) of a particle are taken to be its position at some time t ■ l 0 , say, i.e. 
(o>,e) - (*,y,*)| ( -i. 

The two methods are equivalent, except that the Lagrangian form gives more 
information: it tells where eaeh fluid came from originally, this can considerably facili- 
tate calculations of flows when two or more fluids, with different equations of state, are 
present. However, the Eulerian method has been the popular one and generally 
regarded as more advantageous than the Lagrangian for the following two reasons. 
Firstly, it directly leads to results showing explicitly the spatial distributions of flow 
quantities, such as the pressure field and the stress field, etc., this is desirable in some 
practical applications, e.g. flow over an airplane. Secondly, and more importantly, for 
steady flow the number of independent variables in the Eulerian method is immediately 
reduced by one - the time variable drops out. This is not easily done in the conven- 
tional Lagrangian method for which four independent variables are apparently still 
needed. This simplicity of the Eulerian method for steady flow over that of Lagrange 
explains why the former is used by the great majority of fluid dynamicists. 

The situation is quite different for unsteady flow for which both methods require 
four independent variables and there seems no obvious advantage of the Eulerian 
method over the Lagrangian. In fact, while there is no distinction between steady and 
unsteady flow in the Lagrangian concept, the Eulerian method treats the steady flow as 
a degenerated case (of unsteady flow). As such, any Eulerian approach that is specially 
designed for steady problems may not readily be generalized to unsteady problems. 
Indeed, it has been demonstrated by Van Dommden and Shen (1080) for unsteady 
boundary layer separation and by Hui and Tobak (1081) and Ilui and Van Roesscl 
(1081) for unsteady high Mach number flow that the Lagrangian method can be more 
advantageous than the Eulerian in treating unsteady flow problems. 

The purpose of this paper is to develop a new aud general Lagrangian method 
which uses as Lagrangian variables three material functions and a Lagrangian time r 
distinct from the Eulerian time t. It will be shown that while the new Lagrangian 


method includes the conventional one ns a special case, it requires ouly three indepen- 
dent variables in describing a three-dimensional steady flow: the Lagrangian time and 
two material functions which for steady flow are stream functions. This then places the 
Lagrangian method on the same footing as the Euierian one, The advantage of the new 
Lagrangian formulation over the conventional one for unsteady flow will also be dis- 
cussed. 


2. THE MATERIAL FUNCTIONS 

2.1 Governing Equations In Euierian Form 

We start from Euierian formulation to derive our new Lagrangian formulation. 
Consider a typical problem of fluid flow resulting from a rigid body moving through it. 
Let the motion of the body be prescribed by the velocity v e (l), of its centre of mass and 
its angular velocity ft(t), and let the fluid be compressible, viscous, homogeneous, and 
isotropic. In a body-fixed coordinate system, the equations expressing conservation of 
mass, momentum and energy are the continuity equation, the Navier-Stokes equation, 
and the energy equation respectively: 

-J- + V(/>v) - 0, (1) 

+ (v • V)v + -Vp - F + ^V(V'v) + 

y*.(~~ + v -vrj + pVv - V(kVT)+ (p + }/rXV v) a , 
with the equation of state 

P “ PUT, (4) 

where p is pressure, p density, T temperature, p the shear viscosity, /» the bulk viscos- 
ity, k the termal conductivity, c f and e v the specific heats of the fluid, and R m e f — e ¥ . 
F in Equation (2) is the sum of the body force per unit mass G and the inertia force 
arising from the motion of the body, 

F « G - + 20 x v + O x (O x x) + ~-x x| (5) 

Let (*(f “ 1 , 2 , 3 ) be body-fixed curviliuear coordinates such that *» 0, say, coincides 
with the body surface. Then equations ( 1 ) to ( 3 ) may be written, with the use of 


( 2 ) 

(3) 



summation convention, 
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F + •*% + | V' V + 7*"7jr ’ * + 

4. iL-j» f *V . r , JsL 4. r* _ r * iHl 

P f 1 *V * *l k ~ *9? 

* (ri,» + ri,r,7 - it,rjt).'J , 

"'•(if * •’'^1 + " ‘ + (* + f")**' 


<«> 


( 7 ) 

( 8 ) 


where the g* 3 's are the components of the inverse to the metric, g is the determinant of 
the metric, the rjj’s are the Christoffel symbols of the second kind, and 

#- w- ^ * (i (V7 »') W 

In the special case of inviscid flow of a perfect gas the equations of motion may be 
obtained from (0) and (7) by putting ft ** /* *» 0, and by replacing (8) by 



(10) 


whore 7 «* e f /e, . 


2.2 The Material Functions 

The key to the development of the new Lagrangian formulation lies in the handling 
of the continuity equation. 

It has long been known that for steady two-dimensional or axisym metric flows a 
stream function may be used to eliminate the continuity equation and reduce the 
number of unknown functions by one. it was Lagrange (1781) who first proposed the 
stream function for the two-dimensional motion of an incompressible fluid and sixty-one 
years later Stokes (1842) used a stream function for axisymmetric flows of an 
incompressible fluid in precisely the same manner. Since then the use of stream func- 
tions has become a standard technique for solving two-dimensional or axisymmetric 
flows whether compressible or incompressible, The use of two stream functions for 
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three-dimensional flows, while less well known, has been considered by a number of peo- 
ple, among (hem Yih (1057). 

For unsteady flows the presence of the term in the contiuuity equation (1) 

however, precludes the use of stream functions. But, facilitated by a change of nota- 
tion, it is possible to introduce a set of new functions to be called material funetiona. 

Thus letting ■ I, and noting that the metric is time independent, we may 
re-write the continuity equation as 


o, ini 

where o # ■ 1 and upper case lating indices range from 0 to 3. We now introduce three 
functions »k , (( 0 ,( , ,( 3 *(*)» • “ 1,2,3, related to the components of velocity as follows 


V7 /-•>'« (/ "0,1,2,31 (12) 

where X is an arbitrary function of the t|»"s. Substitution into (11) shows that the con- 
tiuuity equation is identically satisfied. Taking the material derivative of i|»' yields 


D** 

Dt 


+ v i**L m v 'i*L 

dt dt> H* 

VJ/ { * } h ! w K *t L 


- o. 


(i - 1,2,3) 


(13) 


This shows that the i^’s are constant following the fluid and so the name material 
functions is appropriate. The existence of these functions was known by Yih (1057) but 
they have been little utilized since then. 

The mnterial functions introduced above for unsteady flow are a generalization of 
the stream functions for steady flow. In the latter case the time variable t(« (°) does 
not appear, consequently only two material functions ik'UVf 3 ,^ 3 ), »*■ 1,2, are needed to 
eliminate the continuity equation, and (12) reduce to 

VJ yo* - e’^M* 1 ,* 3 )-^ (f “ 1,2,3) (14) 

which show that the material functions in steady flow are just the well-known stream 
functions. 
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Sine* ^ » const, (i • 1,2,3) following n fluid particle, n particle can, in turn, be 
identified by n act of the const ante i|»*, In other words, the material functions t|t* can be 
.and will be in this paper, used as Lagrangian variables. We also note that the material 
functioiM, analogues to the stream functions in steady flow, are not uniquely determined 
by (12) but, rather, there are freedoms in choosing them. The conventional choice of 
the Lagrangian variable (a,t,c) as the particle's position at some time t 0 >* just * pur* 
titular choice of the material functions. In this regard, the conventional Lagrangian 
formulation is a special case of the present one. 

It should also be noted that although Yih (1057) pointed out that the arbitrary 
function may be taken to be uuity without loss of generality, it is not always 
desirable to do so in practical applications as illustrated by an example in section 0. 

The body surface, which must be a material surface, can be described by 4>' ■» 0, 

say. 


3. THE LAGRANGIAN TIME 

We now further introduce a function 4^4° »$\£ 5 .{*) such that 

<*.**.) X(*>V) J‘ 


(15) 


When {!{! m 0,1, 2, 3) are regarded as functions of t| i j [J ■ 0,1, 2, 3), it can easily be 
shown that 


9 I* I 

“ « ■» v 1 . 


( 10 ) 


i,e. 


and 


it 

«s|»° 


* u' (» ■ W). 


(17) 


(18) 


Also, the material derivative reduces to 


D m a 

Dt " a* 0 ’ 


( 10 ) 


which show, that ia a measure of time following a fluid particle or Lagra ngian time. 


♦ 
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For thi* reason, +* will be replaced by r whenever convenient. From now on we shall 
use i |/(/ m 0,1, 2, 3) as our Lagrangian variables. 

The Lagrangian time r differs from the Eulerian time I by an arbitrary function of 

i ■ r + *(<►', <Ji 9 ,ik s ), (20) 


as easily seen from Equation (17). This arbitrary function » may be taken to be sero, 
thereby identifying the two times in the unsteady flow case. However, one can also 
take advantage of this new freedom in applications. For instance, in the case of inviscid 
supersonic flow with shock wave past a body, one can choose a(4»*) to be such that the 
unknown shock wave position is given by r •• 0. This amounts to Fixing the sero time 
of a particle at the moment when it crosses the shock. As shown by Van Roessel and 
Hui (1080) this renders the free boundary problem of flow with shock wave in the 
Eulerian {('} space to a fixed boundary one iu the Lagrangian {$') space, the latter 
being easier to solve. 


In the case of steady flow, while the Euleriau time 1 docs not appear, the Lagran* 
gian time r(»i{»®) ‘till does and is defined by 


W, <W) ' 


( 21 ) 


where the two stream functions t|f ( and t|» 9 ore defined by (14). Together with these 
stream functions, (r,4» , ,t|r) form a set of three independent Lagrangian variables for 
describing a thrce*dimcnsional steady flow. Here r plays the dual role of being the time 
variable and of distinguishing fluid particles. 


The Lagrangian time is thus seen to be more fundamental than the Eulerian time 
when both steady and unsteady flow are considered. 


4. LAGRANGIAN FORMULATION 

Making the transformation of independent variables from {{'} to {«k'} in (6) to (8), 
one obtains the following equations of motion in Lagraugian form 

-jjr + *"*%■*£■ * - r + 



r 






dr 


dr 


+ (w,|rj,) + r;,r/j - rj,r;i)-^} , (.' - 1.2,3) 

* 

/f.~ + pi - ni|i'|| + (»■ + ~p)t\ 

Hi ,<', <*,<’) . M*‘, ♦*,♦*) 

*(r.+‘. *»,*>) V7, 


( 22 ) 

(23) 

(24) 


I - r + .(t'W) 


(25) 


where 



and the operator JV/ (/ ■ 0,1 ,2,3) is 


Ni - 



VJp [ ag y ai K dj L d dj* dj K dj L d 
X e, ’ Wi ldi>‘ d^ 2 dll'* ~dr ' d Iff 3 dtk* dr dll' 1 . 

*■ 

+ *t' dt K d( L d _ 

dll' 1 dr d^ 1 dll' 2 ” dr dll' 1 dll' 2 dll' 3 j 


(20) 


(27) 


Equation (22) to (25) are six equations for six unknowns p,p,{ f (/ ** 0,1, 2, 3). Of 
course, Equation (25) can be used to eliminate I and (24) to eliminate p, resulting in 
four equations for four unknowns. 

Two-dimensional and one-dimensional flow equations in Lagrangian form arc given 
in the Appendix. 


6* SPECIAL CASES 


6.1 Invbcld Flow 

In the case of inviscid flow of a perfect gas Equations (22) are simplified by putting 
#*“/**» 0, while (23) is replaced by 



which it readily integrated. 


5.2 Steady Plow 

For steady flow, the two stream functions f l rnd tjf 1 and X(t|i , ,i|r) are denned by 
(14), whereas the Lagrangian time r is defined by (21). The governing equations for 
steady flow in Lagrangian form can be derived by a method (Van Roessel, 1085) similar 
to that leading to Equations (22) • (25) for uustendy flow. Alternatively, and whieb 
amounts to the same, the steady flow equations can be obtained as a special case of the 
unsteady flow equations by contracting both Eulerian time f and the third material 
functions t|i* to the Lagrangian time f. Thus Equation (24) reduces to 



while Equation (25) becomes redundant. On the other hand, Equations (22) and (23) 
remniu unchanged in form except that the operator N, becomes 


* X d$* dr a<|»* dr di^ 1 dr d4*’ di|i 2 / 


(30) 


5.3 Cartesian Coordinates 

In cartesian coordinates ? ■ (x.y.c), ■ 0 g** *» A' 1 , g ■ 1, Equation (22) and 

(23) reduce to 

■0- + - r + (( - 1,2,3) (31) 


»r. . 
'‘•77 + rt 


N,|lr.V,|r|| + (jr + 


(32) 


whereas (24) to (27) are unchanged in form except that } B I, 
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6. AN EXAMPLE! GERSTNER WAVE 


Consider two-dimensional, iuviscid, incompressible, flow of wafer under gravity # . 
Use cartels coordinates ay, with s boritonlal and y vertically upward. Let f - r and 
♦' - o, «y* ■ fi. The governing equations are 


JkiFl . MsJi 

•W) P * 


* 9 * . JL A(m1 

* X *W) 


o, 


i!l + 1 

•»» x *W) 


-y. 


It is known that this system admits an exact analytic solution representing travel* 
ling waves (Gerstner, 1802) as follows 

x « o-oc*^sin(krt - «r/, 

y ■ ft + oc 4 ^cos(ko - wr), 

i 

f m “pfl 8 «V*'-l), 

X ■ />( I - a*kV**), 

where a and k are arbitrary constants, and w 8 “ yk. 

It is noted that X » X(/f) in this case. If X ■ const was imposed, the same 
Gerstner wave solution would be extremely difficult to obtain. 


7. CONCLUSIONS 

t 

Starting from the Gulerian formulation, a new and general Lagrangian formulation 
for fluid flow is developed. This new formulation uses as Lagrangian variables three 
material functions and a Lagrangian time distinct from the Gulerian time. 

Whilst the Gulerian time is the same for all fluid particles, the Lagrangian time 
differs for different particles. It thus plays a dual role of being the time variable and 
simultaneously of distinguishing fluid particles. With the use of Lagrangian time, 
steady flow can be described by three independent variables - the Lagrangian time and 
two stream functions, thus placing the Lagrangian formulation for steady flow on the 
same footing as Gulerian. This is iu direct contrast with the conventional Lagrangian 
method which apparently still- requires four independent variables for describing a 



steady flow - the Eulerian time plus three variables to identify n particle in three* 
dimensional space. 

For unsteady flow the new formulation incl*;-' the conventional formulation as a 
special case when the Lagrangian time is identified with the Eulerian time and when the 
material functions are taken to be the particle position at some given time. However, 
the distinction between Lagrangian and Eulerian time may prove useful in practical 
applications. 

Finally, the present formulation allows simultaneous presence of both Lagrangian 
and Eulerian variables, thus facilitating the change from the former to the latter to give 
results showing tbs explicit spatial distribution of flow quantities. An example of this is 
given by Van Roessel and llui ( 1080 ). 
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APPENDIX; LOWER-DIMENSIONAL FORMULATION 

For two«dimensiona) unsteady flow, Equations (22) • (25) become 

y + y»*'Af,|p| - f° + £ ^-i‘‘ , N 0 |l| 

+ - r;,AU^| 

+ (w,|r;,| + rr,r‘, - r; k ri,) <«- u) 

/K.-J7 + p* - ^=-Af.|V7t,“'w,ini + [?+ }p]#>, 


(Al) 

(as) 


»(•■?.(•) . A(»S» 8 ) 
«(r»4»'»4»*) y/gp 

I « r + e(«J>\i|r), ' 


where 


N , - fi£ + i£ + _»_1 

x a* 8 ar * a^ 8 ar ai|i* ar at|>* a<j* 8 /’ 

O'- 8 0,1,2) 


(A3) 

(A4) 

(A5) 


I 



(AO) 


Similarly, for one*dimcnsioual unsteady flow, rj| ■ •£— 
■ (|», iV| ■ N, and the governing equations reduce to 




J, f" - F, (< 




/ \2 

+ i wip| f + 


4* 


JL 

M 


.«■# + *l££l). 


(A7) 


ar , 
' e ’77 + ' # 


^=-N(^;V|ril + (/i + J71)* 8 , 


(A8) 
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- >(») 
#(T * W G P 

(A9) 

t - T + •(♦) 

(A10) 

N . 0 f 3t 3 3t 

A (at 3* " 3i|» 3rJ * 

(All) 

0 ■ ~ N[/g |£] 

(A12) 
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